Reduced coupling between offline neural replay events and default mode network activation in schizophrenia

Schizophrenia is characterized by an abnormal resting state and default mode network brain activity. However, despite intense study, the mechanisms linking default mode network dynamics to neural computation remain elusive. During rest, sequential hippocampal reactivations, known as 'replay', are played out within default mode network activation windows, highlighting a potential role of replay-default mode network coupling in memory consolidation and model-based mental simulation. Here, we test a hypothesis of reduced replay-default mode network coupling in schizophrenia, using magnetoencephalography and a non-spatial sequence learning task designed to elicit off-task (i.e. resting state) neural replay. Participants with a diagnosis of schizophrenia (n = 28, mean age 28.2 years, range 20-40, 6 females, 13 not taking antipsychotic medication) and non-clinical control participants (n = 29, mean age 28.1 years, range 18-45, 6 females, matched at group level for age, intelligence quotient, gender, years in education and working memory) underwent a magnetoencephalography scan both during task completion and during a post-task resting state session. We used neural decoding to infer the time course of default mode network activation (time-delay embedding hidden Markov model) and spontaneous neural replay (temporally delayed linear modelling) in resting state magnetoencephalography data. Using multiple regression, we then quantified the extent to which default mode network activation was uniquely predicted by replay events that recapitulated the learned task sequences (i.e. 'task-relevant' replay-default mode network coupling). In control participants, replay-default mode network coupling was augmented following sequence learning, an augmentation that was specific for replay of task-relevant (i.e. learned) state transitions. This task-relevant replay-default mode network coupling effect was significantly reduced in schizophrenia (t(52) = 3.93, P = 0.018). Task-relevant replay-default mode network coupling predicted memory maintenance of learned sequences (ρ(52) = 0.31, P = 0.02). Importantly, reduced task-relevant replay-default mode network coupling in schizophrenia was not explained by differential replay or altered default mode network dynamics between groups nor by reference to antipsychotic exposure. Finally, task-relevant replay-default mode network coupling during rest correlated with stimulus-evoked default mode network modulation as measured in a separate task session. In the context of a proposed functional role of replay-default mode network coupling, our findings shed light on the functional significance of default mode network abnormalities in schizophrenia and provide for a consilience between task-based and resting state default mode network findings in this disorder

Post-stroke upper limb recovery is correlated with dynamic resting-state network connectivity

Motor recovery is still limited for people with stroke especially those with greater functional impairments. In order to improve outcome, we need to understand more about the mechanisms underpinning recovery. Task-unbiased, blood-flow independent post-stroke neural activity can be acquired from resting brain electrophysiological recordings, and offers substantial promise to investigate physiological mechanisms, but behaviourally-relevant features of resting-state sensorimotor network dynamics have not yet been identified. Thirty-seven people with subcortical ischemic stroke and unilateral hand paresis of any degree were longitudinally evaluated at 3 weeks (early subacute) and 12 weeks (late subacute) after stroke. Resting-state magnetoencephalography and clinical scores of motor function were recorded and compared with matched controls. Magnetoencephalography data were decomposed using a data-driven Hidden Markov Model into 10 time-varying resting-state networks. People with stroke showed statistically significantly improved Action Research Arm Test and Fugl-Meyer upper extremity scores between 3 weeks and 12 weeks after stroke (both p < 0.001). Hidden Markov Model analysis revealed a primarily alpha-band ipsilesional resting-state sensorimotor network which had a significantly increased life-time (the average time elapsed between entering and exiting the network) and fractional occupancy (the occupied percentage among all networks) at 3 weeks after stroke when compared to controls. The life-time of the ipsilesional resting-state sensorimotor network positively correlated with concurrent motor scores in people with stroke who had not fully recovered. Specifically, this relationship was observed only in ipsilesional rather in contralesional sensorimotor network, default mode network or visual network. The ipsilesional sensorimotor network metrics were not significantly different from controls at 12 weeks after stroke. The increased recruitment of alpha-band ipsilesional resting-state sensorimotor network at subacute stroke served as functionally correlated biomarkers exclusively in people with stroke with not fully recovered hand paresis, plausibly reflecting functional motor recovery processes

Bayesian optimization of large-scale biophysical networks

The relationship between structure and function in the human brain is well established, but not yet well characterised. Large-scale biophysical models allow us to investigate this relationship, by leveraging structural information (e.g. derived from diffusion tractography) in order to couple dynamical models of local neuronal activity into networks of interacting regions distributed across the cortex. In practice however, these models are difficult to parametrise, and their simulation is often delicate and computationally expensive. This undermines the experimental aspect of scientific modelling, and stands in the way of comparing different parametrisations, network architectures, or models in general, with confidence. Here, we advocate the use of Bayesian optimisation for assessing the capabilities of biophysical network models, given a set of desired properties (e.g. band-specific functional connectivity); and in turn the use of this assessment as a principled basis for incremental modelling and model comparison. We adapt an optimisation method designed to cope with costly, high-dimensional, non-convex problems, and demonstrate its use and effectiveness. Using five parameters controlling key aspects of our model, we find that this method is able to converge to regions of high functional similarity with real MEG data, with very few samples given the number of parameters, without getting stuck in local extrema, and while building and exploiting a map of uncertainty defined smoothly across the parameter space. We compare the results obtained using different methods of structural connectivity estimation from diffusion tractography, and find that one method leads to better simulations

osl-dynamics, a toolbox for modeling fast dynamic brain activity

Neural activity contains rich spatiotemporal structure that corresponds to cognition. This includes oscillatory bursting and dynamic activity that span across networks of brain regions, all of which can occur on timescales of tens of milliseconds. While these processes can be accessed through brain recordings and imaging, modeling them presents methodological challenges due to their fast and transient nature. Furthermore, the exact timing and duration of interesting cognitive events are often a priori unknown. Here, we present the OHBA Software Library Dynamics Toolbox (osl-dynamics), a Python-based package that can identify and describe recurrent dynamics in functional neuroimaging data on timescales as fast as tens of milliseconds. At its core are machine learning generative models that are able to adapt to the data and learn the timing, as well as the spatial and spectral characteristics, of brain activity with few assumptions. osl-dynamics incorporates state-of-the-art approaches that can be, and have been, used to elucidate brain dynamics in a wide range of data types, including magneto/electroencephalography, functional magnetic resonance imaging, invasive local field potential recordings, and electrocorticography. It also provides novel summary measures of brain dynamics that can be used to inform our understanding of cognition, behavior, and disease. We hope osl-dynamics will further our understanding of brain function, through its ability to enhance the modeling of fast dynamic processes

Integrating cross-frequency and within band functional networks in resting-state MEG: A multi-layer network approach

Neuronal oscillations exist across a broad frequency spectrum, and are thought to provide a mechanism of interaction between spatially separated brain regions. Since ongoing mental activity necessitates the simultaneous formation of multiple networks, it seems likely that the brain employs interactions within multiple frequency bands, as well as cross-frequency coupling, to support such networks. Here, we propose a multi-layer network framework that elucidates this pan-spectral picture of network interactions. Our network consists of multiple layers (frequency-band specific networks) that influence each other via inter-layer (cross-frequency) coupling. Applying this model to MEG resting-state data and using envelope correlations as connectivity metric, we demonstrate strong dependency between within layer structure and inter-layer coupling, indicating that networks obtained in different frequency bands do not act as independent entities. More specifically, our results suggest that frequency band specific networks are characterised by a common structure seen across all layers, superimposed by layer specific connectivity, and inter-layer coupling is most strongly associated with this common mode. Finally, using a biophysical model, we demonstrate that there are two regimes of multi-layer network behaviour; one in which different layers are independent and a second in which they operate highly dependent. Results suggest that the healthy human brain operates at the transition point between these regimes, allowing for integration and segregation between layers. Overall, our observations show that a complete picture of global brain network connectivity requires integration of connectivity patterns across the full frequency spectrum

How Sensitive Are Conventional MEG Functional Connectivity Metrics With Sliding Windows to Detect Genuine Fluctuations in Dynamic Functional Connectivity?

Despite advances in the field of dynamic connectivity, fixed sliding window approaches for the detection of fluctuations in functional connectivity are still widely used. The use of conventional connectivity metrics in conjunction with a fixed sliding window comes with the arbitrariness of the chosen window lengths. In this paper we use multivariate autoregressive and neural mass models with a-priori defined ground truths to systematically analyse the sensitivity of conventional metrics in combination with different window lengths to detect genuine fluctuations in connectivity for various underlying state durations. Metrics of interest are the coherence, imaginary coherence, phase lag index, phase locking value and the amplitude envelope correlation. We performed analysis for two nodes and at the network level. We demonstrate that these metrics show indeed higher variability for genuine temporal fluctuations in connectivity compared to a static connectivity state superimposed by noise. Overall, the error of the connectivity estimates themselves decreases for longer state durations (order of seconds), while correlations of the connectivity fluctuations with the ground truth was higher for longer state durations. In general, metrics, in combination with a sliding window, perform poorly for very short state durations. Increasing the SNR of the system only leads to a moderate improvement. In addition, at the network level, only longer window widths were sufficient to detect plausible resting state networks that matched the underlying ground truth, especially for the phase locking value, amplitude envelope correlation and coherence. The length of these longer window widths did not necessarily correspond to the underlying state durations. For short window widths resting state network connectivity patterns could not be retrieved. We conclude that fixed sliding window approaches for connectivity can detect modulations of connectivity, but mostly if the underlying dynamics operate on moderate to slow timescales. In practice, this can be a drawback, as state durations can vary significantly in empirical data

Post-stimulus beta responses are modulated by task duration

Modulation of beta-band neural oscillations during and following movement is a robust marker of brain function. In particular, the post-movement beta rebound (PMBR), which occurs on movement cessation, has been related to inhibition and connectivity in the healthy brain, and is perturbed in disease. However, to realise the potential of the PMBR as a biomarker, its modulation by task parameters must be characterised and its functional role determined. Here, we used MEG to image brain electrophysiology during and after a grip-force task, with the aim to characterise how task duration, in the form of an isometric contraction, modulates beta responses. Fourteen participants exerted a 30% maximum voluntary grip-force for 2, 5 and 10 s. Our results showed that the amplitude of the PMBR is modulated by task duration, with increasing duration significantly reducing PMBR amplitude and increasing its time-to-peak. No variation in the amplitude of the movement related beta decrease (MRBD) with task duration was observed. To gain insight into what may underlie these trial-averaged results, we used a Hidden Markov Model to identify the individual trial dynamics of a brain network encompassing bilateral sensorimotor areas. The rapidly evolving dynamics of this network demonstrated similar variation with task parameters to the ‘classical’ rebound, and we show that the modulation of the PMBR can be well-described in terms of increased frequency of beta events on a millisecond timescale rather than modulation of beta amplitude during this time period. Our results add to the emerging picture that, in the case of a carefully controlled paradigm, beta modulation can be systematically controlled by task parameters and such control can reveal new information as to the processes that generate the average beta timecourse. These findings will support design of clinically relevant paradigms and analysis pipelines in future use of the PMBR as a marker of neuropathology

The psychological correlates of distinct neural states occurring during wakeful rest

When unoccupied by an explicit external task, humans engage in a wide range of different types of self-generated thinking. These are often unrelated to the immediate environment and have unique psychological features. Although contemporary perspectives on ongoing thought recognise the heterogeneity of these self-generated states, we lack both a clear understanding of how to classify the specific states, and how they can be mapped empirically. In the current study, we capitalise on advances in machine learning that allow continuous neural data to be divided into a set of distinct temporally re-occurring patterns, or states. We applied this technique to a large set of resting state data in which we also acquired retrospective descriptions of the participants' experiences during the scan. We found that two of the identified states were predictive of patterns of thinking at rest. One state highlighted a pattern of neural activity commonly seen during demanding tasks, and the time individuals spent in this state was associated with descriptions of experience focused on problem solving in the future. A second state was associated with patterns of activity that are commonly seen under less demanding conditions, and the time spent in it was linked to reports of intrusive thoughts about the past. Finally, we found that these two neural states tended to fall at either end of a neural hierarchy that is thought to reflect the brain's response to cognitive demands. Together, these results demonstrate that approaches which take advantage of time-varying changes in neural function can play an important role in understanding the repertoire of self-generated states. Moreover, they establish that important features of self-generated ongoing experience are related to variation along a similar vein to those seen when the brain responds to cognitive task demands